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INTRODUCTION 


Laboratories  using  a  varied  collection  of  ventilators  report  contradictory 
results  in  animal#  and  nan  regarding  ideal  oscillatory  frequencies  and  stroke 
volumes  to  be  used  in  high-frequency  ventilation  (HFV)  (Bohn  et  al.,  1980; 
Rossing  et  al.,  1981;  Schmid,  Knopp,  and  Rehder,  1981;  Slutsky  et  al.,  1980). 

We  have  examined  the  relationships  between  frequency,  acoustic  power,  and  gas 
mixing  in  a  two-compartment  mechanical  model  of  a  ventilator-lung  system.  This 
model  allows  the  monitoring  of  pressures  and  gas  composition  at  the  site  where 
gas  mixing  occurs. 

The  above  studies  had  two  principle  aims.  The  first  was  to  determine 
whether  power  delivery  to  the  site  of  mixing  controls  mixing  efficiency.  Many 
combinations  of  frequency  and  amplitude  could  in  principle  lead  to  the  same 
power  when  measured  in  the  lung.  The  second  was  to  evaluate  the  effectiveness 
of  random  noise  (excitation  over  a  broad  band  of  frequencies)  in  promoting  gas 
mixing.  If  mixing  were  aonotonically  related  to  power,  and  noise  could  be  used 
for  excitation,  then  the  achievement  of  satisfactory  settings  on  ventilators 
might  be  simplified. 

METHODS 

The  model  (Fig.  1)  was  composed  of  two  compartments,  an  acrylic  speaker 
housing  (compartment  Cl)  25  cm  in  diameter,  with  a  volume  of  5  liters,  and  a 
glass  bottle  (compartment  C2)  15  cm  in  diameter,  yith  a  volume  of  4  liters. 

The  two  compartments  were  joined  by  a  cylinder  19  cm  long  and  3  cm  in  diameter. 

C2  was  always  filled  with  air  at  the  beginning  of  each  experiment,  while 
Cl  was  filled  with  helium,  oxygen,  or  sulphur  hexafluoride.  The  compartments 
were  separated  during  flushing,  and  rejoined  at  the  beginning  of  each 
experimental  run.  Once  the  compartments  were  connected,  gas  was  sampled 
continuously  from  the  distal  portion  of  C2  and  analyzed  for  content  with  a 
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Fig.  1.  Experimental  setup.  For  description  see  text. 

Abbreviations:  Cl  -  compartment  1,  with  the 

speaker;  C2  -  compartment  2;  M.S.  -  mass 
spectrometer;  DVOM  -  digital  voltmeter;  AMP  - 
pressure  transducer  amplifiers;  FFT  -  Fast 
Fourier  Transform  Analyzer;  POP  11/34  -  DEC 
computer;  Ps  -  Valldyne  pressure  monitoring 
for  compartment  Cl. 
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Medspec  II  mass  spectrometer  (Chemetron  Corp.,  St.  Louis,  MO),  or  by  a  Med 
Science  SOS  Nitralyzer  (Med  Science  Electronics,  St.  Louis,  MO)  when  0^  wee 
used  as  the  washout  gas.  The  nitrogen  analyzer  had  a  sampling  rate  of 
3  ml /min,  the  mass  spectrometer  100  ml/min.  The  nitrogen  concentration  in  C2 
was  recorded  on  a  Nicolet  Explorer  II  digital  oscilloscope  (Nicolet  Scientific 
Corp.,  Northvale,  NJ)  and  later  transferred  to  DEC  Nll/70  and  / 34  computers  for 
model  and  statistical  analyses.  The  speaker  was  activated  15  sec  after  the  two 
compartments  were  connected.  Each  experiment  was  3.4  min  in  duration. 

The  speaker  was  driven  by  either  a  sinusoidal  source  (Hewlett-Packard  3311 
Function  Generator)  or  by  pseudo-random  binary  noise  generated  by  a  Nicolet  660 
Fast  Fourier  Transform  Analyzer  (Nicolet  Scientific,  Northvale,  NJ)  and 
filtered  by  a  Krohn-hite  3320  digital  filter  (Krohn-hite,  Cambridge,  MA) . 
Exciting  voltages  were  amplified  by  a  power  amplifier  and,  in  the  case  of 
noise,  by  a  Gould  DC  amplifier  (Gould  Inc.,  Cleveland,  OH)  as  well.  The 
driving  voltages  were  then  directed  to  the  coils  of  a  40-W,  12-ia  loudspeaker 
(Oaktron,  Monroe,  WI) .  The  speaker  diaphragm  was  coated  with  a  thin  layer  of 
latex  to  reduce  gas  leaks;  this  coating  raised  the  speaker  resonant  frequency 
from  35  to  161  Hz.  Elimination  of  all  leaks  could  not  be  accomplished  without 
destroying  the  function  of  the  speaker. 

A  digital  voltmeter  read  RMS  voltage  across  the  speaker  inputs.  Pressures 
were  measured  by  two  5F  Millar  PC-350  Mikro-Tip  Catheter  Pressure  Transducers 
excited  by  TC-100  Millar  Transducer  Control  Units  (Millar  Inc.,  Houston,  TX) . 
These  units  were  in  turn  driven  by  a  transducer  amplifier  (Gould  Inc., 
Cleveland,  OH).  While  one  Millar  remained  fixed  in  the  distal  portion  of  C2, 
the  other  could  be  used  to  monitor  pressures  throughout  the  model.  The  fixed 
Millar  was  adjacent  to  the  lumen  of  the  gas  sampling  catheter  in  the  mid-distal 
region  of  C2.  A  Statham  PM131TC  pressure  transducer  (Gould  Inc.,  Cleveland, 
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OH),  with  diaphragm  exposed,  was  used  in  Cl.  The  transfer  function  describing 
the  response  of  these  transducers  was  flat  within  10X  over  the  range  from  30  to 
R0  Hz. 


The  voltages  applied  to  the  speaker  and  the  resulting  pressure  signals  in 

Cl  and  C2  were  used  by  the  FFT  analyzer  to  calculate  and  display  power  spectra 

(Bracewell,  1978),  power  expressed  as  a  function  of  frequency.  Average 

acoustical  power  (P  )  could  be  determined  from  the  following  relationship 
av 

(Symon,  1964): 


P 

av 


(Ap^)av 
(p  B)-5 


<1> 


where  AP  =  the  amplitude  of  a  pressure  wave,  p  ■  gas  density  and  B  »  the  bulk 
modulus  for  the  medium  through  which  the  wave  is  propagating.  Pgv  has  units  of 
energy  per  unit  area  per  sec  traveling  in  the  direction  n.  With  all  else 
constant,  acoustic  power  calculations  made  by  the  FFT  relied  solely  on  AP  as 
measured  by  the  appropriate  transducer. 

Power  spectra  were  averaged  from  16  samples,  each  sample  consisting  of  4 
sec  of  pressure  recording.  Total  power  (P^,QT),  representing  a  summation  of  the 
power  at  each  of  1046  points  spread  equally  across  a  frequency  range  from  0  to 
100  Hz,  was  used  to  compare  random  and  sinusoidal  inputs.  Transfer  functions 
(Hab)  for  the  physical  system  between  two  transducers  were  calculated  by  the 
analyzer  from  the  power  spectra  recorded  in  both  compartments.  These  functions 
are  expressed  in  the  following  equation: 

Hab  "  Gah/Gaa  <2> 

where  Gab  is  the  average  cross-spectrum  between  the  input  a(t)  and  output  b(t). 
The  cross-spectrum,  rather  than  the  power  spectrum,  of  b(t)  must  he  used  in  the 
numerator  to  insure  that  the  derived  transfer  function  represents  a  linear 
system;  that  is,  the  measured  output  is  due  only  to  the  known  input  and  not 
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some  ocher  unobserved  variable.  This  distinction  is  most  important  when 
transfer  functions  for  the  respiratory  system,  a  notftlly  alinear  system,  are 
obtained.  Caa  is  the  power  spectrum  of  a(t)  (Bracewell,  1978).  Pressure 
signals  were  also  processed  by  a  PDP-11/34  computer  for  time  averaging. 

A  Marquardt  least-squares  regression  (Marquardt,  1963)  was  used  to  fit 
the  gas  mixing  data  to  a  mathematical  model  consisting  of  equations  for  a 
two-compartment,  open  system  (Appendix). 

RESULTS 

Fig.  2  shows  the  change  in  nitrogen  concentration  in  C2,  when  Cl  contained 
02  and  mixing  was  Induced  by  operating  the  speaker.  In  this  example,  Cl 
initially  contained  only  oxygen,  and  the  speaker  was  driven  at  51  Hz  at  80%  of 
the  maximum  power  tolerated  by  the  speaker.  The  foam  suspension  for  the 
speaker  diaphragm  was  permeable,  accounting  for  the  rise  in  N2  near  the  end  of 
some  washouts.  Our  mathematical  model  included  an  expression  for  the  loss  (K^) 
of  tracer  (02  or  SF6)  from  Cl. 

For  a  set  power  and  driving  frequency,  nitrogen  mixing  curves  were 
super lmposable.  When  the  speaker  remained  off,  mixing  induced  by  both 
diffusion  and  the  convective  pull  of  the  gas  analyzers  was  relatively  small 
(Fig.  3).  When  fixed  RMS  voltages  were  applied  to  the  speaker  at  a  succession 
of  frequencies,  an  optimum  frequency  was  found.  A  smoke  plume  within  C2 
revealed  extreme  turbulence  when  the  speaker  was  driven  at  the  optimal 
frequency. 

The  coefficient  of  transport  (T^)  from  Cl  to  C2  was  estimated  by  a 
least-squares  fit  to  the  mixing  data  (Appendix).  These  coefficients  are 
plotted  in  Fig.  4  against  various  frequencies  and  powers  in  Cl.  For  a  given 
coefficient  the  required  input  power  varied  widely  with  changing  frequencies. 

At  70  Hz,  for  example,  voltages  that  would  have  destroyed  the  speaker  would 
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A  nitrogen  mixing  curve  with  gas  sampling  from 
C2,  using  sinusoidal  excitation  at  50  Hz.  0^ 
was  the  diluting  gas.  The  data  and  mathematical 
fit  for  a  two-compartment,  open  system  are 
superimposable. 
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Fig.  3.  Mixing  curves  in  an  0  -N-  environment  with  a 

fixed  driving  voltage  of  10  V  RMS.  An  optimal 
frequency  of  SO  Hz  is  evident. 
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have  been  required  to  achieve  appreciable  mixing. 

Pronounced  frequency  dependencies  were  due  to  the  shepe  of  the  system 
transfer  functions  (Fig.  5).  These  functions  were  obtained  by  applying  random 
noise  over  the  frequency  range  from  1  to  100  Hs,  and  comparing  the  resultant 
pressures  in  Cl  and  C2.  Increasing  gas  density  caused  a  progressive  reduction 
in  resonant  frequency  from  63  to  35  Hs.  The  pressure  wevee  peasing  through  the 
connector  were  so  effective  neer  the  raaonant  frequency  that  mixing  stimulated 
by  other  frequencies  seemed  slight  in  cosq>arison. 

Overall  system  transfer  functions  relating  speaker  driving  voltages  to  the 
pressures  in  C2  were  almost  idee deal  to  those  illustrated  in  Fig.  5, 
indicating  that  the  impedance  of  the  speaker  bed  little  to  do  with  the  overall 
system  properties  for  the  frequency  range  examined.  Measurement  of  the 
impedance  across  C2  was  attested  by  placing  the  movable  Millar  pressure 
transducer  at  the  proximal  end  of  C2  and  comparing  its  output  to  that  from  the 
fixed,  distal  Millar.  This  impedance  was  too  low  to  be  measured;  the  transfer 
function  approximated  1  for  the  frequency  range  1  -  100  Hz. 

When  coefficients  for  mixing  of  nitrogen  and  oxygen  were  plotted  against 
power  measured  in  C2  (Fig.  6),  we  found  that  the  most  effective  frequencies 
were  the  higher  ones.  For  a  given  coefficient  of  mixing,  however,  the  power 
required  to  achieve  those  mixing  rates  varied  less  with  frequency  than  in  the 
case  of  powers  measured  in  Cl. 

When  total  powers  (PT0T)  in  C2  were  matched  (Fig.  7),  and  02  was  initially 
in  Cl,  the  initial  rate  of  loss  of  from  C2  was  similar  for  both  random  noise 
(40  to  70  Hz  band-pass  filtered)  and  sinusoidal  excitation  at  49  Hz.  (Fig.  8). 
The  same  was  true  when  SF6  and  air  were  mixed  by  applying  matched  powers  for 
noise  and  sinusoidal  excitation  (Fig.  9).  Nevertheless,  sinusoidal  forcing 
resulted  in  a  delay  in  initiation  of  mixing  with  SF6  that  was  2  sec  longer  than 
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Fig.  5.  The  influence  of  gas  density  on  resonant 

frequencies  is  illustrated  by  the  transfer 
function  (TF  =*  Hab  of  Eq  .  1)  relating  C2  to 
Cl.  The  gases  used  in  Cl  were,  from  left  to 
right,  SF6 ,  02,  and  He. 
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The  power  spectra  measured  in  C2  under  both 
sinusoidal  and  noise  excitation  are  compared 
The  noise  is  bandpass  filtered  from  40  to 
70  Hz.  The  total  powers  in  C2  (P  Q_,  the 
summation  of  power  across  frequencies  from 
0  to  100  Hz)  were  made  identical  for  both 
forms  of  excitation. 
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Fig. 


8.  Filtered  noise  and  sinusoidal  excitation  produce 
washouts  that  are  similar  under  conditions  of 
matched  total  power  in  C2. 
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when  random  noise  was  used,  [lateacy^^^  ■  15.7  +  0.5  s  (n*3);  latency>^n>  “ 
17.5  +  0.4  •  (n  ■  4) j  P  <  0.005].  The  X  nitrogen  curve  did  not  rise  during 
these  tests,  presumably  due  to  the  poor  dlffuelblllty  of  SF6  through  leaks  In 
Cl. 

DISCUSSION 

Since  the  Impedance  across  C2  wee  Immeasurably  smell,  the  power  dissipated 
in  C2  must  have  been  smell.  Nevertheless,  the  result  of  that  power 
dissipation,  convective  miming,  wee  dearly  dependent  on  the  power  delivered  to 
C2. 

While  ve  would  have  liked  to  extend  the  power  ranges  measured  in  C2  for 
frequencies  other  than  40  end  50  Hr  (Fig.  6),  this  was  not  possible  due  to  the 
restrictions  on  power  input  to  the  speaker  (40  V)  and  the  narrow  band-pass 
characteristics  of  the  connector.  Nevertheless,  for  low  rates  of  mixing, 
matching  power  at  th~  site  of  mixing  removes  much  of  the  frequency  dependency. 
For  example,  to  maintain  a  given  rate  of  gas  mixing  (T£  -  50  sec-1)  while 
shifting  from  50  to  30  Hz  (Fig.  4),  power  in  Cl  had  to  increase  more  then  20 
times.  In  contrast,  the  power  required  la  C2  changed  no  more  than  twofold  for 
the  same  change  in  frequency  (Fig.  6). 

When  power  was  matched  in  C2,  mixing  generally  was  more  rapid  at  higher 
frequencies,  just  as  would  be  predicted  if  we  assume  peak  flows  followed  peak 
pressures  (Slutsky,  1980).  We  did  note,  however,  that  for  a  given  power  In  C2, 
40  Hz  resulted  in  faster  mixing  than  did  50  Hz.  We  heve  no  explanation  for 
this  observation. 

The  reported  optimal  frequencies  for  HFV  (Bohn  et  el.,  1980;  Rossing  et 
al.,  1981)  may  simply  result  from  equipment  and  anatomical  impedances.  A 
similar  suggestion  has  been  made  by  Schmid,  Kaopp,  end  Rehder,  1981.  Since 
impedance  is  a  function  of  ventilator  design,  tubing  geometry,  and  gas  density, 


15 


as  well  as  anatomy,  a  universal,  ideal  oscillatory  frequency  should  not  be 
expected. 

There  remains  a  need  to  determine  the  effectiveness  of  HFV  in  some  rapid, 
convenient  manner.  This  work  suggests  that  the  measurement  of  airway 
oscillatory  pressures  as  close  to  the  site  of  gas  exchange  as  possible  may 
yield  the  needed  Information.  In  animals,  a  piexoelectric  catheter  transducer 
can  be  inserted  down  an  endotracheal  tube,  thus  allowing  high  frequency 
measurements  in  the  lover  trachea.  It  may  be  that  accelerometer  measurements 
on  the  chest  wall  may  eventually  yield  information  on  the  power  delivered  to 
alveoli,  once  due  account  is  made  of  chest  wall  impedances  and  energy 
transmission  from  the  larger  airways.  Nitrogen  washout  curves  may  be  a  further 
means  of  quickly  assessing  the  effectiveness  of  HFV. 

Random  noise  exhibits  several  useful  characteristics.  If  the  resonant 
properties  of  a  ventilatory  circuit  are  such  that  driving  frequency  is 
critical,  and  that  frequency  is  apt  to  change  for  one  reason  or  another,  than 
random  noise  provides  Increased  latitude  in  setting  ventilatory  parameters. 

If  different  paths  through  the  pulmonary  tree  have  different  resonant 
properties,  the  use  of  a  broad  band  of  frequencies  may  be  f*r  more  effective  in 
promoting  gas  exchange  than  any  single  frequency  at  the  same  input  power. 

Another  property  of  random  noise  is  illustrated  in  Fig.  10.  Records  of 
random  and  sinusoidal  pressure  wave  forms  are  shown  for  the  condition  of 
matched  total  power  in  C2.  The  sinusoid  had  peak  to  peak  pressure  excursions 
of  15  cm  H20.  The  random  noise  had  more  pressure  swings  for  the  same  time 
interval  (more  than  three  times),  but  the  mean  pressure  excursion  was  only  5  cm 
HjO.  Only  22  of  those  excursions  had  pressures  above  12  cm  HjO,  and  only  15% 
had  pressures  greater  than  10  cm  HjO.  Indeed,  it  is  the  nature  of  sinusoids 
that  their  most  probable  amplitudes  are  their  extremes,  whereas  the  most 
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probable  amplitudes  for  noise  are  close  to  the  median,  in  this  case,  zero 
(Bracewell,  197 8).  If  one  wishes  to  minimize  the  mechanical  affacta 
of  oscillations  on  airways  by  minimizing  the  magnitude  of  pressure  changes, 
random  noise  would  seem  to  have  a  clear  advantage  over  sinusoidal  forcing. 

Power,  or  something  akin  to  it,  can  be  monitored  quite  simply.  An  audio 
range  voltmeter  provides  RMS  voltage  that  la  proportional  to  the  square  root  of 
power  for  sinusoidal  signals.  For  random  noise  inputs,  an  AC  voltmeter  was  not 
useful;  instead  ve  used  the  time  average  of  squared  voltages  from  the  pressure 
transducers.  The  time  average  correlated  well  with  power  found  by  averaging  in 
the  frequency  domain,  and  involved  simpler  equipment  and  computations. 

SUMMARY 

In  our  model  we  saw  nothing  to  refute  the  concepts  expressed  by  Bohn  et 
al.,  1980,  relating  enhancement  of  diffusibillty  to  increasing  frequency 
monotonlcally  (Fig.  7).  In  our  model,  however,  the  first  requirement  was  for 
power  delivery  to  the  site  of  mixing.  Power  delivery  was  in  turn  affected  by 
impedances  between  the  oscillator  and  the  mixing  site.  When  there  is  neither  a 
detailed  knowledge  of  system  Impedances  nor  measurements  of  power  at  the  site 
of  mixing,  inferences  about  frequency  dependence  of  high-frequency  ventilation 
must  be  suspect. 

Lastly,  questions  of  optimal  frequency  may  be  sidestepped  by  the  use  of 
noise  for  power  gas  mixing.  Certain  theoretical  and  practical  advantages 
accrue  from  this  mode  of  mixing. 
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Appendix  1 


Riggs,  1963,  provides  equations  derived  from  Skinner  et  al.,  1959, 

describing  the  appearance  of  a  tracer  (02  or  SF6  in  our  case)  in  the  second  of 

two  compartments  (C2),  after  introducing  the  tracer  to  the  first  compartment 

(Cl) .  The  tracer  is  transported  to  C2  at  a  rate  determined  by  the  transport 

coefficient  (T^).  (Transport  from  C2  to  Cl  was  assumed  to  be  governed  by  the 

same  coefficient.)  The  derived  values  for  T  would  have  been  somewhat  lower 

c 

and  related  in  a  complicated  fashion  with  power  If  we  had  not  taken  into 
account  the  leakage  of  oxygen  from  Cl.  The  rise  of  nitrogen  concentration  in 
C2  during  the  latter  part  of  experimental  runs  with  oxygen  was  influenced  by 
both  (the  leak  in  Cl)  and  T^.  Least-squares  regression  was  used  to  solve 
the  following  equation  for  both  and  . 

z  n2  -  (i  -  fc2  t)  x  100 

where  Fc2  t  is  the  fraction  of  tracer  appearing  with  time  in  C2. 

Fc2  t  is  expressed  as: 

Fc2.t  ■!_£.(  [-1/2(2  Tc  +  Kx  -  G) t]  [-1/2(2  T.  +  -  G)t] 

G  (e  -e 

where  G  «  -\Jh2  +  4<V2 

When  is  zero,  the  above  equation  reduces  to  a  single  exponential. 


